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Electric Fields Produced in Cubic Crystals by Point Defects 
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Charged point defects in crystals polarize the surrounding ions. These induced dipoles contribute 
to the electric field in the crystal. In this paper, lattice sums for the calculation of this contribution 
are given for the NaCl, CsCl, CaF 2 , and ZnS structures. Various positions for the defect are chosen, 
and the field evaluated near lattice sites in the vicinity of the defect, with radial (with respect to the 
defect) displacements of ± 20 percent of the cation-anion distance allowed away from each lattice site. 
The lattice sums are expressed in power series, including terms up to cubic in the displacements, and 
the coefficients tabulated. 
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Electric Fields From Defects in Cubic Crystals 

Impurities and other point defects in ionic crystals 
polarize the ions around them. The electric field 
produced at a neighboring ion is the sum of the field 
due to the charge on the defect plus the fields due to 
dipoles (electronic and displacement) induced upon 
all the other surrounding ions. In the dielectric 
continuum model of Mott and Littleton, 2 the induced 
dipole on the ith ion is 



If the defect is taken as the origin, and with rj\rj\ = t> j 
and the vectors expressed as 



r = a(/ix + ky + /z) 



(3) 



where x, y, and z are the coordinate axes, eq (2) can 
be written as 
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q\ The quantities Sjc and Sja are lattice sums over the 
cation and anion sites respectively, and have the form 



where r, is the vector from the defect to ith ion, A<? is 
the difference between the charge on the defect and 
that normally associated with the site occupied by 
the defect, rjj is the volume of crystal containing one 
molecule MaXb, and 



M+ = 



OL± 



K-\ 



aa + -\-ba- 4>ttK 



with a ± the polarizability of the ] ca lon [ an( j /£ t h e 

I anion J 
static dielectric constant of the crystal. 

If the vector from the ith to the 7th ion be designated 
r ih then the field at the ;th ion due to the induced di- 
poles on all the other ions is 
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with 
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p M =hi+q+i1' 



P^hjhi + kjki+ljh 
P«=Af + A?+lf. 

The sum does not include the origin nor the evalua- 
tion point j, but does include in principle all other lat- 
tice sites. The evaluation point j need not be at a 
lattice site; by taking nonintegral values of hj, kj, and 
/j, the sums appropriate to other positions may be 
found. 

The lattice sums have been evaluated by computer 
for several cubic structures, with the defect occupying 
various lattice sites and with the evaluation of the 
field also at or near several different lattice sites. The 
evaluation point was allowed to move away from the 
lattice site a distance and direction specified in each 
case by a single parameter X. This parameter is so 



defined that if r is the radius vector from the defect 
at the origin to the lattice site, (1 + AX) r is the radius 
vector of the evaluation point, with A chosen to make 
the displacement of the evaluation point from the 
lattice site approximately a/10 when X=\. The 
sums were evaluated near a given lattice site with 
X equal to — 2, — 1, 0, 1, and 2. From these five 
values for Sj, the coefficients in the interpolation 
formula 



Table la. Sets of coefficients for use in eq (5) 



Si = S?+XS'.+ X 2 S" + X 3 S" ' 



(5) 



were found by making the series agree exactly with 
the calculated sums atJf= — 1, 0, and 1, and fit ap- 
proximately, by least squares, a.\X = — 2 and 2. The 
values of Sj at \X\ =2 are accurate to about 1 percent. 
Convergence in the calculated fields was assumed to 
have been achieved when the total contribution from 
a given set of ions, with the same numerical values of 
h, k, /, was less than the total from all ions included 
divided by 10 5 . 



Set 


s° 


S' 


S" 


S'" 


A 


-5.8952 


0.6078 


0.2513 


-0.01102 


B 


-1.2153 


- .5796 


- .0599 


.00930 


C 


-2.3713 


.2846 


-.0107 


- .00091 


D 


-0.3880 


- .0538 


.0012 


.00011 


E 


-1.9650 


-.1714 


.0593 


.00755 


F 


-2.2874 


.5183 


.0530 


.00171 


G 


-1.3537 


.1015 


.0020 


- .00066 


H 


-0.7302 


- .0820 


.0050 


- .00019 


I 


- .2467 


- .0305 


.00006 


.00016 


J 


-8.2036 


1.2307 


-.0110 


-.01745 


K 


-0.4697 


-0.1865 


- .0336 


- .00298 


L 


-7.1106 


.0282 


.1914 


- .00173 


M 


-11.5685 


2.3025 


- .0697 


- .03626 


N 


-0.0614 


-0.3048 


- .0935 


- .00829 





- 2.3529 


- .2252 


.0605 


.00766 


P 


- 2.0496 


.8115 


-.1610 


.02028 





- 2.6663 


-.1731 


.0557 


.00189 


R 


-0.9769 


-.1125 


.0050 


.00003 


S 


-2.3570 


.3536 


- .0255 


.00042 


T 


-0.5922 


- .0733 


- .0068 


.00072 


U 


-5.8465 


.8771 


.0145 


-.01787 


V 


0.1222 


-.1132 


- .0404 


- .00370 



Table lb. Electric fields in cubic crystals from dipoles induced by point defects, using eqs 4 and 5 and the coefficients in table la 



Structure 


Distance denoted 
by "a" (eq 3) 


Defect site 


Evaluation point 


Anions 


Cations 




Inside 


Same 
Shell 


Outside 


Inside 


Same 
Shell 


Outside 


NaCl 


Cation-anion 


Cation [anion at (0, 0, 1)] 


(I+O.IXMV2, V2, Va) 
(1+0.08XM0. 0. 1) 
(1+0.05XM0, 1. 1) 


F 


c 


A 
D 
H 




G 


B 

E 

I 


CsCl 


Cation-cation 


Cation [anion at (V2, V2, V2)] 


(I+O.IXXV2, V2, V2) 
U+0.08XM0, 0, 1) 
U+0.05X)(0, 1. 1) 


M 
P 


J 


K 

N 
Q 


F 


c 

G 


L 

R 


CaF, 


Anion-anion 


Cation [anion at (V2, V2, V2)] 


(lH-0.1X)(Va, V 2) V2) 
(1 + 0.8X)(0, 0, 1) 
(1+0.05XH0, 1, 1) 


M 
P 


J 


K 

N 
Q 




G 


B 
E 

I 






Empty cube center 
[Anion at (V2, V2, V2) | 


(I+0.IXXV2, V2, Va) 
(1+0.08XM0. 0. 1) 
(1+0.05XX0, 1. 1) 


M 
P 


J 


K 

N 

Q 


F 


c 


A 
D 
H 






Anion [cation at (V2, V2, V2)] 


(l+0.1X)(Va, V2, V2) 
(1+0.08X)(0. 0, 1) 
(l + 0.LX)(Ya, V2,-V2) 




c 


L 


L 


M/2 


S 

u 


T 

N/2 
V 


Zinc-blende 


(Cation-cation)/ V 2 


Anion [cation at (V2, V2, Va)] 


<1+0.1X)(V 2 . V2, Va) 
(1+0.08X)(0, 0, 1) 

(1+0.1XHV2, 1 /2,- 1 / 2 ) 






B 
E 
B 


M/2 


s 

u 


T 

N/2 

V 






Midpoint of unit cell edge 
[cation at (V2, V2, V2)] 
[anion at (1. 0. 0)] 


(l+0.1X)(Va, Va, Va) 
(1+0.08XX0, 0, 1) 
(1+O.lXXVa, Va,-Va) 




c 


A 
D 
A 


M/2 


s 

u 


T 

N/2 

V 



The results, in the form of the coefficients in the 
power series, are given in table 1. Each set of co- 
efficients occurs several times, so for the sake of 
economy these sets are listed together in table la and 
given symbols. In table lb, these symbols are used 
to indicate the set of coefficients appropriate to the 
situation described in the first four columns of the 
table. 

With the defect at the origin, the ions in the crystal 
can be divided into three groups: those on the shell 



that includes the lattice site around which the evalua- 
tion is being performed ("same shell"), all those within 
this shell ("inside"), and all those outside of it ("out- 
side"). In the last six columns, the contributions to 
Sj from anions and cations in these three groups are 
listed separately. Values for S° symbolized by Z), E, 
r, and V have already been given by Mott and Littleton 
(see footnote 2). 

(Paper 71A1-128) 



